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Abstract 

We consider the gluon fusion production cross section of a scalar Higgs boson in 
models where fermion and scalar massive colored particles are present. We report 
analytic expressions for the matrix elements of gg — ► Hg, qq — > Hg, and qg — > Hq 
processes completing the calculation of the NLO QCD corrections in these extended 
scenarios. The formulas are written in a complete general case, allowing a flexible use 
for different theoretical models. Applications of our results to two different models are 
presented: i) a model in which the SM Higgs sector is augmented by a weak doublet 
scalar in the SU(N C ) adjoint representation, ii) The MSSM, in the limit of neglecting 
the gluino contribution to the cross section. 
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1 Introduction 



The Higgs mechanism for the electroweak symmetry breaking is the still untested part of 
the Standard Model (SM). The search for the Higgs boson is one of the most important 
goals of the present experimental program at the Tevatron and, in the near future, at the 
Large Hadron Collider (LHC) at CERN. 

Electroweak (EW) precision physics data and the direct search limit from LEP constrain 
the possible values of the Higgs mass in the SM quite strongly, with a solid indication that 
a SM Higgs boson should be lighter than 200 GeV. However, in many extentions of the 
SM the above bound does not apply because of the presence of new physics (NP) that can 
affect the EW fit allowing a higher value for the Higgs mass. Theoretical arguments based 
on perturbative unitarity, triviality and fine-tuning indicate that the crucial mass range to 
be investigated is up to the TeV scale. The search of the Higgs boson at the LHC has, 
therefore, to be supported by an accurate theoretical knowledge of the production cross 
sections, the decay modes, and the important background processes in this range (for a 
general review see Ref. pQ). 

Due to the gluon luminosity, the main production mechanism for a scalar Higgs boson 
at the LHC is the gluon fusion process (pp — * H + X) [2j. Its cross section in the SM is, 
all over the range of interesting values of the Higgs mass, one order of magnitude bigger 
than that of the other main production mechanisms, the Vector Boson Fusion (VBF) [3], 
and the production in association with heavy quarks and vector bosons [3]. 

Differently from the other production mechanisms, the gluon fusion is a process that 
starts at 0{a 2 s G lJ ), i.e. at the one-loop level. In fact, the Higgs boson does not couple 
to gluons directly, but only via a loop of colored particles. Thus, the gluon fusion process 
is the Higgs production mechanism where NP can play the most relevant role changing 
significantly the value of the production cross section. 

The first predictions for the gluon fusion Higgs boson production cross section in the 
SM dates in the late seventies [2]. More than fifteen years later, in 1991, the calculation 
of the NLO QCD corrections was completed in the infinite top mass, m t , limit [5] and, 
successively, retaining the full dependence on the mass of the heavy fermion that runs in 
the loops [B]. The total effect of the NLO QCD corrections is the increase of the LO cross 
section by a factor 1.5—1.7, giving a residual renormalization/factorization scale dependence 
of about 30%. The unexpected size of the NLO QCD radiative corrections, that seem to 
spoil the validity of perturbation theory, motivated, at the beginning of 2000, the calcula- 
tion of the NNLO QCD corrections, performed in the infinite mt limit [Sj. The calculation 
shows a good convergence of the perturbative series. The NNLO corrections are sizable, 
but, nevertheless, smaller that the NLO ones. Moreover, the QCD bands of variation of 
the renormalization/factorization scale overlap with the ones of the NLO calculation. The 
NNLO corrections enhance the cross section of an additional 15% (of the NLO results). 
Moreover, they improve the stability against renormalization/factorization scale variations. 
Furthermore, the effect due to the resummation of soft-gluon radiation at the NNLL accu- 
racy has been evaluated in Ref. [5]. Besides an additional enhancement of the cross section 
of the order of some percents (up to 6%), the effect mainly lies in a strong reduction of the 
scale dependence. The remaining theoretical uncertainty, due to higher-order QCD correc- 
tions, has been estimated to be smaller than 10%. This estimate was confirmed recently 
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by the NNNLO calculation of Ref. [TDJ. 

Because of the high accuracy reached in the evaluation of the QCD corrections, also the 
EW NLO corrections to gluon fusion were recently taken into account. In Ref. [11], they 
were evaluated in the infinite m t limit, giving a correction of less than 1%. In Ref. [12], 
the contributions coming from Feynman diagrams with a closed loop of light fermions 
were calculated in a closed analytic form, expressing the formulas in terms of generalized 
harmonic polylogarithms (GHPLs) [13]. It turned out that they are sizeable. In particular, 
in the intermediate Higgs mass range, from 114 GeV up the the 2mw threshold, these 
corrections increase the LO cross section by an amount of 4 to 9%. For m# > 2mw, 
they change sign and reduce the LO cross section; however, in this region the light-fermion 
corrections are quite small, reaching at most a -2%. In Ref. [2], also the remaining EW 
corrections due to the top quark were calculated as a Taylor expansion in m\j{^m\ v ). 
This result is valid in the m# < 2 my/ range in which the corrections due to the top quark 
have opposite sign with respect to the light-fermion contribution. However, the former are 
smaller in size, reaching at most a 15% of the latter. 

Several efforts have also been devoted to the calculation of radiative corrections (mainly 
QCD) to less inclusive quantities, like the transverse momentum (q T ) distribution [T5J [TBI 
rm [T8| [T9] and the rapidity distribution [201 HH [19] . All these results have been implemented 
in two Monte Carlos that calculate fully-differential distributions at the NNLO [22[ [23]. 

The Higgs boson gluon fusion production cross section was also extensively studied in 
the Minimal Supersymmetric Standard Model (MSSM). The effects of the squarks on the 
production cross section for the neutral CP-even h and H Higgs bosons were considered in 
Ref. [23], including the NLO QCD corrections evaluated in the heavy squark mass limit. 
This approximation has been relaxed in Ref. [25] where the full dependence on the squarks 
masses has been retained. The complete MSSM NLO QCD corrections to the Higgs bosons 
production, in the heavy SUSY particles limit, has been presented in Ref. [25]. In general, 
these corrections lead to a NLO K factor that differs from the corresponding SM one 
by an amount less than 5%, with the exception of regions where the squark and quark 
contributions interfere negatively giving rise to a MSSM production cross section much 
smaller than the SM one [27]. Not inclusive quantities were also studied in the MSMM [28J. 

As already pointed out, the Higgs boson production via gluon fusion, as well as the 
Higgs decay into two photons, are processes sensitive to any kind of NP. Thus, it would 
be desirable to have predictions for these quantities, at the level of NLO QCD corrections, 
provided in full generality, i.e. not constrained by a particular theoretical model but flexible 
enough to be used for different models, and, if possible, expressed in an analytic form easy 
to be evaluated numerically. 

The aim of this paper is to analyze in a, as much as possible, model independent way 
the contribution of colored scalar particles to the gluon fusion Higgs boson production cross 
section at the NLO level in the QCD corrections. In this spirit, we present here general 
analytic formulas for the NLO corrections to the production cross section of a Higgs boson 
via gluon fusion, a(pp — > H + X). In particular, we provide analytic expressions for the 
NLO QCD corrections to the partonic processes gg — > Hg, qq — > Hg, and qg — > Hq, in the 
two cases in which a fermion or a scalar run in the loops. Together with the analytic results 
of Refs. [291 SO], where the two-loop virtual QCD correction to the gluon fusion process 
where evaluated, the present work completes the NLO calculation of the Higgs production 
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cross section in the presence of colored scalar particles. 

The paper is organized as follows. In Section [2j we provide all the analytic formulas and 
we discuss the validity of the heavy fermion/scalar mass limit. Applications of our results 
to two different models are presented in the following section. In Section [37T1 we consider a 
model, proposed by Manohar and Wise [3T] , in which the Standard Model is supplemented 
by an additional scalar colored isospin doublet. We discuss, at the NLO level, the effects 
of these additional scalar particles on the production cross section of the standard Higgs 
boson and on its decay width H — ► 77. In Section 13.21 we consider the contributions due 
to the squarks in the MSSM. We focus on the NLO QCD corrections (exchange of gluons), 
neglecting the contributions coming from the gluino. Finally, in Section H] we present our 
conclusions. 



2 Higgs Production via Gluon Fusion at NLO 

In this section we present analytic results for the NLO QCD corrections to Higgs boson 
production via the gluon fusion mechanism. Being the Higgs boson neutral under SU(N C ), 
its coupling to the gluons is mediated by a loop of colored particles. To discuss the gluon 
fusion mechanism in a general way we assume as colored particles one fermion and one 
scalar in a generic R1/2, Ro SU(N C ) representation, respectively. The extension to more 
fermions or scalars is trivial. The coupling's strengths of these particles to the Higgs are 
assumed to be: 

A 2 



HFF = g Ai /5 



™>l/2 

2m„ 



HSS = g An. 



(1) 



where g is the SU{2) coupling, m w is the W mass, m 1 / 2 is the fermion mass, A is a generic 
coupling with the dimension of mass and Aj are numerical coefficient^. 

The hadronic cross section for the Higgs production via gluon fusion at center-of-mass 
energy y/s, can be written as: 



a{h 1 + h 2 ^H + X) 



a,b 



X 



dx 1 d,X 2 faM&^A) fb,h 2 ( X 2i A) X 

dz 8 ( z — ) a ab (z) , 
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where t h = m 2 H /s, fj,p is the factorization scale, f a ,hi(x, fi 2 F ), the parton density of the 
colliding hadron hi for the parton of type a, (a = g, q, q) and a a b the cross section for the 
partonic subprocess ab — > H + X at the center-of-mass energy s = x±x 2 s = m 2 H /z. The 
latter can be written as: 

a ab (z) = a {0) zG ab (z), (3) 

where 
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is the Born-level contribution with mo the mass of the scalar particle, 
Gffl = -Ay 1/2 [2 - (1 - Ay 1/2 ) H(0, 0, x 1/2 
g£ l) = iy [l + 2y H(0,0,x )} 



(5) 
(6) 



and T(Ri) is the matrix normalization factor of the Ri representation (T(R) = 1/2 for the 
fundamental representation of SU(N C ), T(R) = N c for the adjoint one). In Eqs. (JMEJ) 
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and, employing the standard notation for the Harmonic Polylogarithms (HPLs), H(0, 0, z) 
labels a HPL of weight 2 that results to 
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Up to NLO terms, we can write 
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We discuss the various NLO contributions. 

i) The gg channel (Eq. (fTTj) ) involves virtual and real corrections. The former, regular- 
ized by the infrared singular part of the cross section gg — > Hg, are displayed in the first 
row of Eq. ( iTTj) where /3q = (11 CU — 4n/ T(Rf) — n s T(R s ))/6 with n/ (n s ) the number of 
active fermion (scalar) flavor in the representation Rf (R s ). The functions Qf 1 ^ containing 
the mass-dependent contribution of the two-loop virtual corrections, can be cast in the 
following form: 
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2 All the analytic continuations are obtained with the replacement — m 2 H 
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where C(Ri) is the Casimir factor of the Ri representation (in particular, for the fundamen- 
tal and the adjoint representations of SU(N C ) we have Cf = (N% — 1)/(2N C ) and Ca = N c , 
respectively). Explicit analytic expressions for g\ 2l) (z = 0, 1/2) given in terms of HPLs 
can be found in Ref. [29]. It should be noticed that Q\ ' R depend upon the choice of 
the renormalized parameters (masses and couplings). In Ref. [29] expressions for Q^ 1 > Cr ' 
with MS or on-shell (OS) parameters are presented. In the case of single heavy fermion, 
i.e. A = 0, become independent of the renormalized mass chosen and goes to the 
well know result —3/2C(Rx/ 2 ) + 5/2 Ca, that can be also obtained via an effective theory 
calculation [5]. The case of a single heavy scalar (Ai/ 2 = 0) is actually more complicated be- 
cause, in general, the coupling of a colored scalar particle to the Higgs boson is not directly 
proportional to the mass of the scalar and therefore different renormalization prescriptions 
for A and mo can be employed. In case an MS prescription is employed both for A and 
mo, the function tends, for large values of m , to the constant value 9/2C(R ) + Ca 
independent on the MS subtraction scale. 

The second row of Eq. ffTTj) contains the non-singular contribution from the real gluon 
emission in the gluon fusion process where 
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is the LO Altarelli-Parisi splitting function. The function lZ gg can be written as 
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where t = —1(1 — z)(l — v), u — —1(1 — z)v, with 
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Furthermore, the functions A 2 and A4 can be cast in the following form: 
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We find 

s 



b 1/2 (s,t,u) = B 1/2 (s,t,u) + -[H(p,0,x 1/2 )-H(p,0,x a 



s s 2 



2 s + u 
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~H 3 (s,u,t) + ~H 3 (t,s,u), (22) 

b (s,t,u) = --B (s,t,u), (23) 

c 1/2 {s, t, u) = C 1/2 (s, t, u) + — *— [H(0, 0, x 1/2 ) - H(0, 0, x s )] + -^—H 3 {s, u, t), (24) 

co(s,t,u) = --C {s,t,u), (25) 



where 



x a = (a = s, t, w) 

v/l-4/a+l 



and 



1 + — ) # (0, 0, a*) + H{0, 0, x s 



.s 



-2 (-^- 2 -1--) [H(0, 0, x.) - H(0, 0, a*)] 
\(s + uy s J 

+l(— + 3]H 3 (s,u,t)-H 3 (t,s,u), (26) 



2 \ s 

Ci(s,t,u) = -2s-2[H(p,0,Xi)-H(0,0,x,))-H 3 (u,s,t). (27) 

In Eqs. ( 122H27]) H(0,x) = ln(x), and the function H 3 is symmetric under the interchange 
of its last two arguments, i.e. H 3 (a,b,c) = H 3 (a,c,b), as can be seen from its integral 
representation: 

r 1 i 

H 3 (a,b,c) = / dx— r — — r- {lnfl — bx(l — x)] + lnfl — cx(l — x)] 

v ' J x{l-x)+a/{bc) 1 L v ;j 

-ln[l-(a + b + c)x(l-x)}}. (28) 

An explicit analytic expression for H 3 (a, b, c) can be found, for instance, in Ref. [T5j. Using 
the notations of [15], we have H 3 (a, b, c) = —W 3 (b, a,c, a + b + cjE 



3 In Ref. [15], in the function W3 the mass of the heavy particle is explicitely written in the integral 
representation. Instead the H3 function has as input parameters "reduced" variables, i.e. Mandelstam 
variables divided by the heavy particle mass. Therefore, in order to have the correct formal expression for 
H3 one has to put in the formulas of Ref. [15] m/ = 1 and consider the variables s, t and u as reduced 
variables. 
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In the case of a single heavy fermion or a single heavy scalar or both fermion and scalar 
heavy K gg -> -11(1 - z) 3 /(Qz). 

ii) The qq — > iff/ annihilation channel, Eq. (|12p . can be written as 



K 



128 z (1 - 2) 
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In the case of a single heavy fermion or a single heavy scalar or both fermion and scalar 
heavy TZ qq -> 32(1 - zfj{TJz). 

iii) Finally we consider the quark-gluon scattering channel, qg — > qH. In Eq. (1131) P gq 
is the LO Altarelli-Parisi splitting function 
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In the case of a single heavy fermion or a single heavy scalar or both fermion and scalar 
heavy K qg -»• 2 */3 - (1 - 2) 2 /^. 
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Figure 1: Exact K-factor normalized to the effective one. The different lines represent a 
single heavy fermion (solid line) and a single heavy scalar whose mass is renormalized in 
the MS (dashed line) or in the on-shell scheme (dashed-dotted line). 



The expressions reported above are, for the fermionic case, in full agreement with the 
known results in the literature [T5| H~6] . They give the exact, i.e. for any value of the 
particle masses, NLO contribution and correspondingly the exact K-factor defined as the 
ratio between the NLO and LO cross sections. It is interesting to compare the value of 
the exact K-factor with the one that can be obtained via an improved effective theory 
calculation. By the latter we mean a result in which the effective NLO cross section is 
obtained by multiplying the exact LO cross section by the 0(a s ) corrections evaluated in 
the heavy particle limit |7J. As discussed above, while for fermions the NLO contribution in 
the limit of heavy mass is independent upon the definition of the renormalized mass used, 
the heavy scalar NLO contribution is actually dependent on the renormalization conditions 
chosen. 

In Fig. [1] we plot the exact NLO K-factor normalized to the effective one as a function 
of m H /rrii, for the case of a single fermion (continuos line) and of a single scalar. For the 
latter we consider two options: i) On-shell condition for mo and MS renormalization of the 
coupling A defined at the MS /i scale /i = m , with A(m ) = m (dash-dotted line), ii) 
MS renormalization both for A and mo with fi = itlq and A(mo) = mo(/i) (dashed line). 
The results presented in the figure have been obtained assuming a hadronic center-of-mass 
energy ^Js = 14 TeV, using the parametrization CTEQ6M [31] to describe the partonic 
content of the proton, and setting the factorization and renormalization scales equal to the 
Higgs boson mass. These choices will be used also in the figures of the following sections. 
From Fig. [T] it appears that in the fermion case the difference between the exact and the 
effective K-factor is at most 10% and that already when m H < Ivrtx^ the difference is below 
1%. Instead, in the scalar case the situation is more complicated. Both cases i) and ii) 
show a spike at the opening of the m H = 2 m threshold. This spike is due to logarithmic 
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and square root singularities present in the two-loop one-particle irreducible (1PI) virtual 
corrections, coming from diagrams in which a scalar self-energy diagram is inserted in a 
one-loop vertex diagram. When the mass of the scalar is renormalized on-shell, the mass- 
counterterm diagrams show a square root singularity that actually cancels the similar one 
coming from the 1PI diagrams, leaving only an unphysical logarithmic singularity (see dash- 
dotted line in Fig. [T]) related to the inadequateness of the standard mass renormalization 
procedure in case of unstable particles. Instead, in the case an MS renormalization for the 
mass of the scalar is employed, the cancellation of the square root singularity between 1PI 
diagrams and mass-counterterm diagrams does not take place anymore. Then both the 
square root and the logarithmic unphysical singularities are left (see dashed line in Fig. [p. 
In the region away from the threshold, the figure shows a good convergence to 1 for light 
Higgs masses. For a heavy Higgs, instead, we note a certain deviation of the effective theory 
from the exact one. Nevertheless, this deviation remains quite limited in size, and it reaches 
at most 9%. 

3 Scalar Particle Effects on the Higgs Production Cross 
Section 

In this section we discuss the effect of colored scalar particles on the Higgs production 
cross section via gluon fusion. We consider two cases: i) a model in which the SM fields 
are supplemented by a weak doublet of colored scalars. ii) The squark contribution in the 
MSSM. 

3.1 The Manohar-Wise Model 

The model proposed by Manohar and Wise (MW) [31] is an extension of the SM that 
includes additional colored scalar fields which transform in the (8, 2)1/2 representation of 
SU (3) x SU (2) x U(l). The choice of these additional scalar fields is dictated by requirement 
of natural suppression of flavor changing neutral currents. The additional colored scalar 
weak doublet 



with a = 1 ... 8 an adjoint color index, contains an electrically charged and two neutral real 
scalars. Denoting the standard (1, 2)1/2 Higgs field by H, the most general potential can 
be written as [3~T] 



where S = S a T a , the trace is over color and SU (2) indices are explicitly shown. The 
ellipses represent tri- and quadrilinear interaction terms of the fields S a not relevant for our 
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Figure 2: Higgs production cross sections as a function of the Higgs mass. The particles 
running in the loop are: only SM fermions (dashed lines) or both SM fermions and MW 
scalars (solid lines). The lower curves in the two cases represent the LO cross section, 
while the upper curves represent the NLO one. The couplings are \i(ms) = 4, ^(ms) = 
1, A 3 (m 5 ) = 1/2. The scalar mass is m s = 750 GeV. 



discussion. The tree-level mass spectrum of the colored octet scalars is found to be: 

2 2 a V 2 

m s+ = m s + A l J , 



™i H = m| + (A 1 + A 2 + 2A 



3 J 



4 

v 2 



ml 0I = m 2 s + (Ai + A 2 - 2A 3 ) — . (39) 



and the coupling to the standard Higgs are: 

A v 2 

HSlS b _ = g -i 5 a 



4 m 



qa qb A X + A 2 + 2A 3 v 2 



o m w 

HSS,S>, = g Al + A ;~ 2A3 . (40) 

o m w 



Eqs. ( 139|40|) show that in this model A/mo ~ v/rris for ms ^> v ensuring the decoupling 
of the colored scalars as their mass increases. 

To analyze the effect on the Higgs production cross section of this octet of colored 
scalars we plot, in Fig. [2, the LO and NLO cross sections, as a function of m H , including 
the scalar contribution and compare them with the SM results. The figure is drawn taking 
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Figure 3: NLO decay width, T 2 l{H —>■ 77) in the SM and in MW model, normalized to the 
LO decay width T\l,sm{H — > 77) in the SM. The dashed line represents the two-loop EW 
SM corrections; the solid line represents the two-loop QCD SM corrections. The dotted 
line represents the full SM two-loop corrections. Finally, the dashed-dotted line represents 
the corrections due to the full two-loop SM+MW decay width. The MW parameters are 
chosen as in Fig. [2] 



T(R Q ) = C(Rq) = 3, m s = 750 GeV, Ai(m s ) = 4, A 2 (m 5 ) = 1, A 3 (m 5 ) = 1/2, (MS 
couplings at the scale fj, = msjj and renormalizing the mass of the scalars on-shell. With 
this set of parameters A 2 /mjj ~ 0.1, thus it acts as a large suppression factor. The figure 
shows that the NLO production cross section in this model is always significantly larger 
that the SM one. In particular, for small values of m H it is almost two times the SM cross 
section. 

The presence of an additional octet of scalar particle affects not only the production 
cross section of the standard Higgs boson but also its decay into two photons, that is a very 
relevant mode for Higgs searches up to m H = 140 GeV. The formulas for the decay width 
H — > 77, including the contribution of colored scalar particles evaluated at the NLO, can 
be found in Ref. [29] . 

In Fig. [3] we plot the correction to the decay width T(H —>■ 77) originating from the 
two-loop (NLO) corrections, with respect to the one-loop (LO) SM prediction, assuming 
the same parameters as in Fig. [2j In the figure, the EW and QCD SM corrections are 
separately shown as well as their sum. As already pointed out in Refs. [321 [33] , the SM 
NLO EW and QCD corrections basically cancel each other. The effect of the charged scalar 
particle, on top of the Standard Model particles, is to reduce the LO SM decay width into 
two photons. At the leading order (one-loop amplitudes), this reduction ranges already 

4 This set of parameters is consistent with the electroweak precision physics constraints [31]. In this 
model it is possible to generate a positive contribution to the p parameter from the colored scalar sector 
allowing for an heavier standard Higgs boson in the electroweak fit. 



11 



between 9 and 6%. Considering the QCD corrections to the scalar contribution (Fig. [3] 
dashed-dotted line), the reduction is further increased by an amount that ranges between 
13% for m H ~ 100 GeV and 6% for m H ~ 160 GeV. The reduction effect clearly depends 
on the mass of the scalar particle and can be much more pronounced for smaller values of 
m s . 

It should be recalled that the relevant quantities at LHC for the Higgs discovery are 
the product of the production cross sections times the branching ratios, so that for m H < 
140 GeV the relevant quantity is a{pp — > H)BR(H — > 77). In this product the reduction 
effect induced in H —>■ 77 by the scalar contribution is actually more than compensated by 
the increase in the gluon fusion production cross section, so that the Higgs boson discovery 
potential at LHC in this model is actually higher than in the SM. 

3.2 The MSSM 

We consider now the contribution of the scalar quarks on the Higgs production cross section 
in the MSSM. The Higgs sector of the MSSM contains two complex (1, 2)1/2, (l)2)_i/ 2 
scalar fields that couple to the down- and up-type fermions separately. After spontaneous 
symmetry breaking the spectrum of the MSSM Higgs sector contains five physical states, 
two CP-even neutral boson, h, H, one CP-odd neutral one, A, and two charged Higgs 
particles H ± . At the lowest order the MSSM Higgs sector can be specified in terms of two 
independent parameters, usually chosen as m^, the mass of the pseudoscalar boson, and 
tan/? = v 2 /vi, the ratio of the vacuum expectation values of the two Higgs fields. 

Stop and sbottom loops can affect significantly the production cross section of the CP- 
even Higgs bosons. Indeed, there can be regions of the SUSY parameter space in which one 
of these squarks can be relatively light and its coupling to the h boson relatively strong. As 
a result this state does not decouple in the gg — > h amplitude and actually its contribution 
interfere with the fermion one. It should be recalled that when going to the NLO level the 
purely squark contribution is only part of the MSSM QCD correction to the production 
cross section. Indeed, at this level, besides diagrams containing quarks or squarks and 
gluons also diagrams with quark, squark and gluino can contribute and a clear separation 
between the two contributions is not possible [25]. However, we are going to consider a 
scenario in which one squark is supposed to be relatively light and therefore to provide 
the bulk of the corrections while gluino diagrams are supposed to give a small contribution 
that we are going to neglected. Also we are not going to take into account the quartic 
self-interaction coupling among squarks. 

The computation of gg — > h/H cross section in the MSSM requires the knowledge of 
the particle mass spectrum of this model. Nowadays, it is available a set of computer codes 
[35] that allow to compute the entire MSSM spectrum starting from a restricted number of 
parameters, that can be assigned at a high scale and then evolved down to the weak scale, 
like for example in a MSUGRA scenario, or directly assigned at the weak scale. At the 
level of NLO corrections it is important to specify exactly the meaning of these parameters. 
For what concerns the entries in the squarks mass matrix, the output of these codes is 
usually expressed in terms of dimensionally reduced DR parameters evaluated at some 
specified \i = Hewsb scale. Consequently the mass eigenvalues obtained from this squarks 
mass matrix, as well as the couplings of the squarks to the neutral Higgs bosons should 
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be intended as DR quantities^]. Among all the various quantities entering in the formulae 
for the gluon fusion cross section at NLO only Q^ l '° R ^ requires an exact specification^. In 
Ref. [29] this quantity is reported in terms of dimensionally regularized MS masses and 
couplings. Thus, to employ the result of Ref. [29] we have first to convert the DR masses 
and couplings obtained as output from any DR code into MS quantities. For what concerns 
the masses one notices that among the various parameters entering in the squark mass 
matrix 

2 _ ( m \ h +m g + m2 z( I q - e g sin 2 6W)cos2/3 m q (A q - fi (cot/?) 27 ^) \ 

q V m q (A q - ji (cot /?) 2/ 9 ) m\ R + m 2 + m| e q sin 2 6 W cos 2/3 / 

the soft SUSY breaking left- and right-handed squark masses, rn~ L , rn~ R , the trilinear squark 
coupling, A q , tan/? as well as the Higgs mass parameter /i at the level of NLO QCD 
corrections are identical in dimensional regularization and dimensional reduction while the 
only parameter that requires a conversion is the quark mass as 



m (DR )=m (MS)__|L CFm9 _ (42) 



In Eq. ( HTj) i| is the third component of the weak isospin, e q the electric charge of the 
quark q, m z the mass of the Z boson and 9w the Weinberg angle. Once an MS squark 
mass matrix has been constructed one can obtain the MS mass eigenvalues and in case 
convert them in the OS results. A similar procedure can be employed to obtain the MS 
couplings of the squarks to the Higgses [361 125], where also in this case the only quantity 
that requires a conversion is the quark mass as in Eq. (j42l . 

Having set the framework for the computation of the gluon fusion production cross sec- 
tion of the neutral MSSM Higgs bosons, we consider the particular region of the parameter 
space, the so-called Higgs gluophobic scenario, in which there is a negative interference be- 
tween the standard fermionic contribution and the one coming from the stop and sbottom 
states [27]. As input parameters for the squark mass matrix at the scale Hewsb = 300 GeV 
in this scenario we chose m~ = m~ = m~ = 350 GeV A t = A b = —600 GeV, u = 300 GeV, 

= 153 GeV, m^S^^) = 2 .3 GeV. 
In Fig. H] we plot the gluon fusion production cross section for the h, H CP-even Higgs 
bosons at LO and NLO for tan (3 = 3. The MS squark mass eigenvalues are found to be: 
m h = 190 GeV, m h = 500 GeV, m~ bl = 350 GeV, m~ h2 = 360 GeV while the rest of the 
MSSM particle spectrum, in particular the masses of the lighter and heavier neutral CP- 
even Higgs bosons, is obtained using the code Suspect with a gluino mass m g = 500 GeV 
and M2 = fi. As can be seen from the figure, when the QCD corrections are taken into 
account the NLO cross section shows an increase comparable to the SM case. Therefore, 
the QCD corrections to the quark and squark contributions are both large and of similar 
size. According to the discussion in Section [2] the NLO curve should contain spikes in 
correspondence of the opening of the 2t 12 , 2 6 12 thresholds. These spikes are actually 
extremely narrow and either are not drawn or are just hinted in the figure. In Fig. [5] the 
same analysis is performed for tan/3 = 30 with a corresponding squark mass spectrum 



5 The codes usually provide also OS masses for the SUSY particles. 

6 For the top and bottom contribution to the gluon fusion cross section we consider always OS masses. 
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Figure 4: Production cross section of a light (heavy) CP-even Higgs boson, in the MSSM, 
with tan/5 = 3 at LO (dashed line) and at NLO (solid line). 
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Figure 5: Production cross section of a light (heavy) CP-even Higgs boson, in the MSSM, 
with tan/3 = 30 at LO (dashed line) and at NLO (solid line). 



14 



m~ tl = 230 GeV, m h = 490 GeV, m~ hx = 320 GeV, m~ b2 = 380 GeV. As can be seen from 
the figure, in this case the NLO corrections are usually percentually smaller than in the 
SM case. Furthermore, in the singular behaviour at the openings of the squarks thresholds 
it is possible to appreciate the change of sign when passing through the thresholds due to 
our choices of MS masses (see Fig. [T]). Our results for the MSSM are in agreement with the 
analysis carried out in Ref. 



4 Conclusions 

In this paper we have presented analytical results for the NLO QCD corrections to the 
Higgs production cross section in gluon fusion. We considered both the contributions due 
to colored fermions and colored scalars running in the loops. The analytic formulas are 
provided in a fully general form and they can be used in computer codes aiming at the 
phenomenological description of different theoretical models. In particular, the results have 
been implemented in a Fortran code which provides a flexible tool to study BSM physics 
effects for a generic model which satisfies SU(N C ) x SU(2)l x U(1)y gauge invariance. 

We have discussed the behaviour of the K factor (the ratio between the NLO and the 
LO cross section), comparing the exact results with the ones in which the NLO corrections 
are calculated in the infinite fermion and/or scalar mass approximation. In particular, we 
recover the well known fact that in the SM the effective theory provides a good approx- 
imation of the NLO corrections in a wide range of values of the Higgs mass. The scalar 
case, however, is more complicated. Actually, a non-physical singularity appears in the 
two-loop virtual corrections at the scalar pair-production threshold (with different shapes 
depending on the renormalization scheme in which coupling and masses are renormalized) . 
This alters, to some extent, the discussion about the behaviour of the K factor near the 
threshold, where the propagator of the scalar field, in principle, should be resummed to all 
orders in perturbation theory. In the region away from threshold, we see a good convergence 
to 1 in the light-Higgs region, while a certain deviation from the effective theory shows up 
for heavy Higgs masses, ranging in any cases within 10%. We notice that, when there is 
a substantial interference between the fermion and scalar contributions the situation could 
be more complicated, giving rise to behaviours that could differ substantially from the ones 
described above. 

In the paper, we have applied our results to the study of the Higgs production cross 
section in two different extensions of the SM: the model proposed by Manohar and Wise, 
in which the SM is supplemented by an extra colored scalar weak doublet in the adjoint 
representation of SU(N C ), and the MSSM in the limit of neglecting the gluino contribution. 
In the MW model, the extra scalars lead to a large enhancement of the Higgs production 
cross section: with the set of parameters considered, we register an enhancement of up to a 
factor of 2 with respect to the SM results. In the same model, the Higgs decay width in two 
photons is decreased by a factor up to 13%. The net effect, considering the combination of 
production and decay, is a large positive correction. In the MSSM, we consider the so-called 
gluophobic scenario in which the destructive interference between squark and quark loops 
reduces significantly the production of the lightest CP-even Higgs boson h. In this situation 
we find that the NLO corrections to the squark contribution are of similar size as those of 
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the quark part in agreement with previous results in the literature [2U [2~5] • 

The study of exclusive observables will be necessary to obtain more realistic phenomeno- 
logical results: the squared matrix elements described in this paper can be easily embedded 
in a Monte Carlo code aiming at such studies. 
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